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Chapter I 
Introduction 

Problems frequently encountered in complex multiple input-ouput systems 
stem from interaction or coupling between the controlled variables. For ex- 
ample, in the lateral control of an aircraft, a change in the angle of bank by 
aileron deflegtion will result in a change in yaw angle provided that there is 
no rudder input to offset the change. To eliminate this interaction between 
the controlled outputs, the theory of M cs controls was developed. 
For a control system designed in accordance with this theory, a desired 
change in one output variable will not cause a change in the other output 
variables. 

The general method for designing such non-interacting control systems 
was first developed by A. S. Boksenbom and R. Hood, Ref. 1, and applied 
by them to the control of a turbo-jet engine. This work was later repeated by 
H. S. Tsien, in Ref. 2. In Ref. 3, H. Freeman extended the theory to obtain 
a desired degree of interaction in an m input and n output system. In Ref. 4, 
Freeman discussed stability and realizability but he did not discuss the practical 
design problems or use the more realistic open loop approach. 

This study will develop the basic criteria for non- interaction of the 
variables in an n input and n output system. In addition, the effect of dis- 
turbance inputs and errors in mathematically describing the physical system 
and synthesizing the controller will be discussed. It will also introduce the 
potential designer of a non-interacting control system to some of the practical 


problems which are often encountered in applying the basic non-interaction 


1 





criteria. As a final step, the theory will be applied to the design of an auto- 
pilot system in order to illustrate the design techniques and the difficulties 


that may be encountered in a practical application. 





Chapter II 
General Theory 
In this chapter, the basic relations between controller and controlled 
elements of a System to effect non- interaction of variables will be developed. 
To achieve this, a specific system configuration will be assumed. However, 
the extension ef this theory to other configurations is feasible and should be 
apparent from the subsequent development. 


Introduction 


A linear physical system with n inputs, 7 , and n outputs, 9 , can be 


== 


simply described in vector notation. In such notation, 7 and d are con- 


sidered as n dimensional vectors with components d; and d, . The coefficients 


of the equations relating 7 to Å in the physical system form an n by n 
matrix, Å , which is composed of elements au Note that the E are oper- 


ators expressable as Em) where p = d/dt. Thus, 


N 


4 = > Uy ہہ‎ 2 (2-1) 


= 
| 


or in vector-matrix notation, 


AS (220)‏ زا 
In a block diagram, the system described above can be represented as:‏ 


Pere ለ በመመ 


In the case where there are m inputs and n outputs and m is greater than n, 


the m-n additional inputs are fed directly through the system to become system 





outputs. This implies that the A matrix is now an m x m matrix with the 
added 0 elements having a value of unity. Such a case is described at 
length in Refs. 1 and 2. 

To achieve non-interaction, it is desired that Y= x where X is an n 
dimensional input vector with components x. In designing a linear system 
to satisfy the above input-output requirements, the foremost problem is to 
use the error vector, € = q-x, to determine ር for non-interaction and at 
the same time ensure that E% 0. One possible system configuration for 
achieving these ends is shown in Fig. ۰٠ 

In order to keep the treatment of the subject sufficiently general, servos, 
transducers and controllers have been included in the overall system shown 
in Fig. l(a). The transducers have been supposed in light of the fact that the 
real outputs, 7 , of a physical system can not be obtained directly. Thus 
the measured outputs, di , are related to the real outputs, Yi , through the 
transducer matrix, T. This matrix is a diagonal matrix with elements t;;. 


It is further supposed that the transducers are selected such that the ሃስ UN 


are very small and hence no attempt will be made to correct di to di : 


The physicalsystem inputs, Di , are generally servo driven and the 7 are 
related to the servo signal inputs, di , through the servo matrix, S. Like 
the T matrix, Sis a diagonal matrix with elements Bu The controller, which 
relates the system input signal vector, ó , to the error signal vector, c 

is represented by an n by n matrix, C. It will be designed such that the 
overall system is non-interacting and has desired response characteristics. 


As a means of simplifying the ensuing mathematical analysis, the system 
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shown in Fig. l(b) shall be used. It should be noted that this is an equivalent 
system to that of Fig. l(a) and that many multiple input-output systems can be 
put into this form. This does not mean that this is the only form possible. 
However, conditions for non-interaction are more complex and system design 
becomes exceedingly difficult, if not impossible, when configurations different 
from Figs. l(aPand 1(b) are used. A simple exception to the above is shown in 
Fig. l(c). Here the diagonal matrix R, which might for example provide deriva- 
tive feedbacks, has been added. Although another loop has been added, the 
system of Fig. 1(c) reduces to that of Fig. 1(b) if A = (T + R)AS and 1 (ፐ+ RIG: 
Analysis of the Control System 


The system under consideration is the simplified but equivalent system 


shown in Fig. 1(b). The vector equation for this system is: 
j= AS = ACe = AC (x-4) (2-3) 
Solving equation (2-3) for the system output vector, 1 , provides: 


4 [Trac] AC X (2-4) 


where I is the identity matrix. 

In order that the inverse matrix [I FAC]. shall exist, it is required 
that the determinant, | I+AC| , is not identically zero over all values of the 
complex variable p. Of course, AC will be a function of the operator p and 
hence زا‎ +۵6 may be equal to zero for a few particular values of the operator 


p considered as a complex variable. 


If the system under consideration is to be non-interacting it is required that: 





qu Dx (2-5) 


where D isa diagonal matrix with elements di. This also implies that 
Ui = di Xj for all i or that a specific output is affected only by a correspondi 
input. = 


From equations (2-4) and (2-5) it follows that: 
-| — 
[E +ac] AC = D (2-6) 


Certainly, one condition on the matrix AC that satisfies the above equation is 


=1 
that AC is 8 diagonal matrix. If AC is diagonal then [r4Ac] AC is given by: 


| (ac), 
I+ (AC) 
2 Ma T voc 
[I+AC] AC = : Is O : (2-7) 
Ou ى‎ ._. O 
(AC); 


0 Qe HAC) ن‎ 
where (ae) are the components of AC. 


Solving equation (2-6) for AC yields: 
اب‎ — 
Ac = [r-5] D (2-8) 


Again, itis obvious that AC is a diagonal matrix if equation (2-8) is to be 
satisfied. 
In summarizing, a necessary and sufficient condition for non-interaction 


is that the matrix AC be diagonal; i.e. AC = D. Furthermore, the elements 





(ac) or d; must be finite and not identically minus one. This assures that: 


۱ e" Chi T : 
Yi 2 1ن‎ E ach _ Xp >= په‎ for L= | ton (2-9) 
HOC) | +d 
where d;; is finite. 
It is now known that AC - D where D is a diagonal matrix different from 


the diagonal matrix D. Solving for C by matrix methods yields: 


C-AD © (2-10) 


For e to exist, the determinant of A must not be identically equal to zero. 
This requirement is also necessary for n inputs to specify n outputs independently. 
Thus, a condition where IA O implies a lack of independence among the 
system equations and the impossibility of completely controlling n outputs Ui 
by n inputs m 


Equation (2-10) can be rewritten as: 
ገገ 
seus -i 
Cu > Lie Ot (2-11) 
Ral 


D is a diagonal matrix with elements di; and therefore the terms in the summation 


are non-zero only for k = j. Thus, 
Cij = كت‎ )2-12( 
Considering the case where j = i, equation (2-11) reduces to Cj; — TG dij om 
di = a (2-13) 
Cu 


Note that equation (2-13) is the open loop transfer function of the overall system. 





7 
Here 0 it is known and 5: is designed so that the closed loop transfer function, 
dii 


/: , provides favorable response characteristics. 
| + مل‎ 





After the c;; have been found by conventional compensation methods, the 
off diagonal elements of the controller matrix, Cjj, are found by combining 


equations (2-12) and (2-13) to obtain: 


| 7 ve 
Jr مه زب‎ 
زد‎ 


If any of the dai =0, then the choice of di is arbitrary. It follows that 
the off diagonal components, e are then determined from equation (2-12). 
Although 0 is arbitrary when A =O, it should be chosen such that good closed 
loop response characteristics are obtained. Thus, it should be chosen fairly 
large in order that outputs approximate inputs with a minimum steady state 
error. This may be seen from the relation, di = du :. Naturally, the d.. 

+9ፈ፪ + ii 

can not be overly large otherwise the 5. must assume impractical proportions 


through equation (2-12). 


In summary, the c,, are determined by the design of the open loop transfer 





function, di = Cu where 01 is known. Thus the c,, are essentially compen- 


ii 





sation so chosen that is stable with favorable response characteristics. 


3 
The off diagonal elements of the controller matrix are then found from equation 
(2-14). Provided no difficulties are encountered in compensating the di the 


system performance should be as desired, namely, non-interaction of control 


variables. 
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Chapter III 
Error Analysis 

There are essentially two basic types of system errors. The first type 
of error is a dynamic error resulting from random disturbances entering the 
System. The second type are those errors introduced into the system through 
the A and C mátrices due to inaccuracies in describing the true physical system 
and/or inaccuracies in the synthesis of the controller elements or transfer 
functions. Errors of the second type could obviously result in some system 
interaction and conceivably stability problems. 
Disturbance Errors 

In analyzing errors of the first type, the disturbances were considered 
as input signals to the A and C matrices as shown in Fig. 2. It is supposed that 
m disturbance signals, VIER are introduced into the C matrix and r disturbance 
Signals, Na , are introduced into the A matrix. These input signals may be 
represented as an m dimensional vector, U , and an r dimensional vector, V : 
As an example of the source of these inputs, the U type of input might be 
considered as N* noise or drift type signals introduced at the inputs of the YO 
“1 components. The ۷ type of signals are introduced into s physical system 
being controlled and hence could be torque loads, gust effects, 5 Regardless 
of their source, it is desirable to minimize the effect of the U and V dis- 
turbances. 

It will be assumed that ٨ and V are related to the system in vector- 


matrix form as follows: 





c J A‏ مسج 


ከ Introauction of Disturbances into the 


Control System 








ፅ= ርፎ ተ ር*ዚ ۱ 
47 AS + AV (3-1) 


In the above equations, C* is an n by m matrix and A* is an n by r matrix. 


Noting that €= X-Y , the system error equation can be found from equations 


(3-1) to be: E 


€ -«[r«AQ] (x-Ac'u- A'v) o 


Specifically, the errors due to the U and V inputs can be expressed as: 
3 

፦ [I+ac] AC*u 
T 

- [I+AC] ۸۷ 


Eu 
Ev 


(3-3) 


(t 


In designing a control system, the designer will have to investigate the 
errors resulting from the U and V disturbances. This will require an investi- 
gation of the above vector equations to determine the error associated with the 
specific input disturbance. This may be illustrated by considering the effect 
of a disturbance input to a single component of the C matrix in a two-dimensional 
System. Assuming the disturbance as an input, U, , to the وه‎ component of 


the C matrix, equation (3-3) becomes: 


Eu = -Ea ۳ الم‎ ٣ 


Or 


Eu, 


- [Ere a P u, 
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As developed previously, for non-interaction, AC is a diagonal matrix which 
was denoted as D. Therefore, 





Eu, = — | t+dy 5 Aula) (41 
O FE (Qa Ciz) 
and 6 27 ۷ Ge Giz U, | 
* | 1+du | 
C, = m 042 Ciz ፌ( 
2 i+d22 


For illustration purposes, 6 will be the only error considered. From 
i 


equations (2-13) and (2-14): 


zi 
dy = EL 1 ns = Qiz Caz = - C,, 
d Gas 1 
xj 
(2 5 T ر)‎ = - = Cu 
Therefore, 
En መ وه‎ 
Ä [+ Cn 317 


Assuming the hypothetical system represented by: 





[yu A 

P p(pt!) 
Am) رر‎ 
pti p*! 


then, 


.! O ዯ ፕ . ብ... لتم د با‎ Mf 
ar daz p(p+1) P' à 
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Now, if the control matrix diagonal elements are designed as: 


Ci = 1 Coz = a 
P 


then, 





Therefore, C 2 | | 
| p (pti ) 


Since |. then the output y, with Xi assumed as zero, will be: 


LL !‏ ع 
TET Sp]‏ 
It is readily apparent that this transfer function presents an undesir able‏ 
situation, since at low frequencies the magnitude of the output will be large‏ 
because of the term .2/p. In addition, the output diverges with time for a‏ 
step input, U, . Thus, in this hypothetical case the designer would be forced‏ 
۱ 

to modify the controller elements in order to provide good syStem response 
and at the same time minimize the disturbance errors, which may be difficult, 
if not impossible. 

One method of reducing the effect of the disturbances which may prove 
effective in a practical design is to introduce a high gain prior to the disturbance 
inputs. Fig. 3(a) illustrates a basic circuit for a scalar system (one of the non- 


interacting subsystems for example) where, with X=Q, the transfer function 


relating the output, 4 , to the disturbance input, UL , is: 





15 


Fig. 3(a ) 





Fig. 3(b) 
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With the introduction of a high gain element A as in Fig. 3(b), the circuit 
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can be reduced to that shown in Fig. 3(c). From this latter figure, with X=0, 


ملد له 
I+ YP A v2‏ 1 


By proper design, {4 can introduce high gain over the desired frequency 
range and consequently reduce the output due to the disturbance inputs. JA 
is, of course, designed to maintain good system characteristics for 7 with 
X. as an input. 
Component Errors 
To analyze the second type of error it will be necessary to refer to 


equation (3-2) which is repeated below: 


€= fr+ad (X-Actu-A*V) a 


The elements of the A, A*, C and C* matrices can be redefined in terms of a 


lowest common denominator as follows: 


|... Qu (p) 
Ee بس‎ 
E Gy (gy 


Cu 


1 
(5 
v, 


dc Cp) (3-6) 


x 


en 2 Ci; (p) 





1% 


where ae ; C : Qe; ; Co are polynomials in the operator p and d, (P) and 
d.(p) are the lowest common denominators of the aij and Ci; components re- 
spectively. 

From the above definitions, it is apparent that the e have the same poles 
as the a by definition. The same is true for E: and E This generally is the 
case and if it is’not, then the nonacceptable poles may be removed and used to 

| 


define a new vector, UL or V . Itis now possible to défine new matrices 


whose elements are polynomials in the operator p. Thus, 

(Qi) 2 

Am ۰ Ax - 

( C ሣ | ) ESE 
E 


Note that A = SA, oe 
G- ۵ 


equation (3-2) may be written as: 


2 


(3-7) 


e 
1 
١ 


CS > غه‎ Hence: 


€= | dade 1 se AC] (dade X- ጳር"ዚ- d, Av) (3-8) 


‘Ne 4 ሥሪ Ya 
In the above equation, the components of Å ; A ; E ; E ; da and de 
are polynomials in the operator p. Hence any poles of equation (3-8) result 


dad. I AC | . If determinant | dude T+ AC | 


from the determinant of 





has no zeroes in the right half p plane, or on the jw axis (imaginary axis), then 
the system is stable for disturbance inputs as well as normal inputs, X; . 
However, if determinant | dod L + AC | and d d. have zeroes which are 
common, and if these zeroes are in the right half p plane then the system will 
be stable for X inputs but could be unstable for U and V inputs. JH dade 


has no zeroes in the right half p plane and the system is designed so that 
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du ۳ 22 
| td, PS |* di; 


inputs, X , U and Y. 


are stable, then the error vector, € , must be stable for all 


When there are small changes in the coefficients of A and C due to inaccur- 
acies in the synthesis of C or inexact knowledge of the coefficients of the equations 


d.d. T T NE will 


of the physical system, the zeroes of the determinant of 





change because Of the coefficient change. Hence, the poles of equation (3-8) 
will be altered slightly for small changes in the coefficients. Thus, if the 
system is designed stable and if d, d. has no zeroes in the right half complex 
plane we are assured that the system remains stable. 

The small changes in the coefficients of A and C will result in interaction. 
However, if the deviations are truly small, then the resulting interaction will 
be negligible. This is due to the fact that the dij or off diagonal elements will 
be small and hence will have little effect on the overall system. Thus, it may 
be seen that interaction due to these slight deviations can be minimized only 
through careful attention to the analytical description of the physical system 


and accurate design and mechanization of the controller. 





Chapter IV 
Practical System Problems 

In the previous chapter, the effect of small deviations, E n aij 
and T were discussed. In this chapter, the problems or consequences 
attendant with gross deviations in ei because of design advantages will be 
investigated. Ig addition, problems arising from impossible or impractical 
compensation requirements, as may be dictated by non-interaction theory, will 
be discussed and means of alleviating the problem pointed out. 
Problems arising from gross deviations in 51 

Suppose that in synthesizing complicated controller transfer functions, 
Cip simplifications are made by neglecting terms in the numerator and 
denominator of : that have "break'' frequencies well above the natural 
frequency of the system. Further suppose that all these discarded terms have 
negative roots. The obvious consequence of such gross deviations from the 
desired Cjj is system interaction at the higher frequencies. In addition, one 
must consider the effect of these departures on the stability of the overall 
System. 

To investigate these ideas more fully, assume a two-dimensional control 


system in which the exact controller transfer functions, c 12 and co ] are 


altered by a factor q to produce new relations Clo and SMS Analytically, 


/ 
Cia = 1 Ciz (4-1) 
a = fs Cer 


(4-2) 
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For the theoretical non-interacting system, the following equations have been 


shown to exist: 


du = Qu Cu + Qz Ca = Y= ay (4-3) 
d, s ue t az عم‎ 2 6, - CS (4-4) 
d, = Qu Ge + hz Ca = O (4-5) 
de = du Cn + Gaz ረ. ZrO (4-6) 


Substituting equations (4-1) and (4-2) into equations (4-3) through (4-6) where 


applicable and manipulating and rearranging where desired, the results are: 


/ = -l 
dy = Au Ci + Quz Ga Car = [Out faa GY | em 


d daz Ca + $ le ۵ TAT da fpa [Ca (4-8) 


a) = Un ٩ Cre +d Cze = D (1-4, ) Gal FO (4-9) 
4 = du Cn + fe lai dg = | (1-42) (u|FO (4-10) 


where the di are elements of a new and non-diagonal matrix, D'. 

The interaction is evident from equations (4-9) and (4-10). The changes 
in the transfer functions d;; due to the factors q; may be determined from 
equations (4-7) and (4-8). 

In investigating the effects of gross deviations on system stability, it is 


best to look at the system equation in vector-matrix notation. It is: 
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Y= E +D] Dx 


where D' is a non-diagonal matrix. 
The stability of the closed loop system is determined from the determinant 
of [I+D | which is the characteristic equation of the overall system. It can be 


shown that: 


|I+D| = (14d, JU +d;2)- 9958, (4-11) 


It follows that if this resulting polynomial has zeroes in the right half complex 
plane the overall system is closed loop unstable. 
One means of determining the stability of the closed loop system without 


solving for the roots of the polynomial is as follows: 


(1+ di |) 1+ ويك‎ ( - diz d2 = 0 


or when rewritten: 
/ 


/ 
(id) (322) ٠ -Sedu — = 0 4-12 
(ltd X141; E 

Stability is assured if none of the above factors have zeroes in the right 
half p plane. 

The open loop transfer functions, di and 922) which are given in 
equations (4-7) and (4-8), should be plotted in Nyquist plots and Nyquist stability 
criteria, N- P- Z applied to determine if there are zeroes of the 1 + dii in 


the right half p plane. 
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Assuming that there are no zeroes of the 1+ di; in the right half p 


plane and that 913 and 05 1 have no poles in the right half p plane, then for 
/ 
a 23, 
(1+9ሀ)(|+ሀ፤,) 


unity for all frequencies, W . However, reference to the Nyquist plot of 


system stability, the absolute value of must be less than 


a typical dii Shown in Fig. 4 provides the following important relationship: 
$ 


| ل + ا‎ | = A where a= 5 (4-13) 
| Mi +1 


It follows from equation (4-13) that: 
/ / / 
diz dzi > لك‎ 
ll v (i tdi roy 


Thus, if (1+ M,) (1 + M,) d'a dh, < 1 for all frequencies, C2 , and (Ir d 


12 1? 


as well as (1 + do o) have no zeroes in the right half p plane, then the overall 


control system is closed loop stable. 


Compensation Problems 


It is conceivable that the co-factors of the 91, of the physical system 


may have zeroes in the right half p plane so that the will have poles in 


| 
ant 
the right half p plane. From the basic open loop non-interaction relation, 
e = Cu , it follows that it may be impossible or impractical to compensate 
il 
AE with Ci; SO that di is stable. In this event, it would be impossible 





Ci HF dii 


to obtain a practical non-interacting control system using the design techniques 


discussed in the theory section. As an example, consider the following simple 





System: = cg 4 
c , OF (5 | 
2 -,Olpt I 2 





pri pt | 
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and the real axis. 


1 
di plane 
Mi aq: The distance from the 
| 1 MT ICO pont to the In- 
| ter section ot the M circle 
| 
| 







| | 
“44 | መሚ | بت‎ — 


EU E کے‎ 
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M circle f 
For MIM, 


/m aginary 5 


ከ ከውር Nyquist Plot of a Typical Open Loop Transfer 


Function. da: 
7 am 





ed a 3 (-.01p+1)(p+1) 


oue T. 46 


and: 


6+ ما٥.-‏ ے كل 
(-Olp+1) (pr)‏ 0 


The locus of = for p = jw is: 
it = plane 





۲۰۱٢٢٣ stability, N= P- Z and in this case P = 1. Therefore, N must be ۴ 1 


if Z is to be zero. However, from the general non-interaction theory, d= ai : 
u 


It is quite obvious that the compensation of xs Such that N= +1 and Z = Û is 


a difficult problem. One means of compensation is to use the prediction oper- 
ator, ef . This will provide the necessary encirclements but is entirely 


impractical as a means of compensation in a control problem. Other than this, 


no compensation could be found. One might consider mechanizing the compensa- 


tion so that it cancelled the objectionable poles and zeroes. However, this will 
not work because the poles and zeroes of the physical system are not known 
exactly and unless we have exact cancellation the poles and zeroes are not 
really cancelled. Thus, it appears that the straight forward design techniques 
developed in the theory can not be used when confronted with such a compen- 


sation problem. 





In some cases the foregoing design problem can be avoided by merely 
interchanging the rows of the ó vector. This implies an interchange of 
controls. If the design problem is removed by such å change, one may continue 
to design the control system in accordance with the general non-interaction 
theory. 

AS an exarfple of the foregoing, consider the same system of the previous 

/ 
example but with the b vector rows interchanged to form a vector ó and a 


corresponding change of the columns of A to form a new matrix A'. Thus, 








AS 
ec | pra +۱ 8 
A = (5) 
-:O1p+I E 
pt! pri 


The new di; are now: 


d= Yo, = Al = ~ Cu C2)Golp +96) 


a pt! 


d,, = 18 5 -= A | - - Ces (2) (— 00 p ዓሪ) 
m pri 
There are no poles in the right half p plane of the above open loop transfer 
functions. Hence, P= 0. In addition, the design of ره‎ and Coo So as to have 
N= Z= 0 in the closed loop system is not anticipated to present any real 
compensation difficulties. Thus, the impossible or impractical compensation 
problem has been avoided. It should be emphasized that interchanging the rows 


of the Ó vector will not alleviate this particular compensation problem in all 
| 


( 
i 
۹ 
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instances. This is obvious from the foregoing example if either Or 


RE 

Aiz 

a, have one or more poles in the right half p plane. Naturally, the presence 
21 


of such poles does not imply that the can not be practically compensated 





ac} 
so that Z = 0 but rather that difficulties E compensation are more probable 
than the instances when P = 0. 
In a more 6omplicated system where there are n control inputs and out- 
puts and n is greater than two, the poles of =r result from the co-factors of 
ul 


the A matrix. These complicated co-factors tend to increase the probability 


| 
that ES will have one or more poles in the right half p plane. Thus, one 


(1 
might expect some difficult and impossible compensation problems when applying 
the general non-interaction theory to complicated physical systems. If an inter- 
change of the controls, (interchange di) does not avoid these problems then a 
different design approach must be taken. 

The theory behind one approach and the effect on system performance is 
most easily described for a simple two-dimensional system. For instance, 
suppose that SET can not be properly compensated by Cil The result is that 
the system can not be designed entirely within the non-interaction theory. How- 
ever, a practical system can be found if the designer will accept some inter- 
action. Furthermore, some of the non-interaction techniques can be used once 
the difficulty has been removed. 


In finding a practical system it is probably best to start with the basic 


expression for di: This is: 


d, = Qi Cy + GizCa, (4-14) 
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In as much as it is impossible or impractical to compensate dii with Ci] 


then c,, may in general be chosen arbitrarily. However, in some instances 651 


11 


may be chosen completely arbitrary and hence the arbitrariness of 511 will be 
restrained to a certain type of transfer function which results in 1 4 di having 


no zeroes in the right half p plane. In addition, certain simplifying assumptions 
| 

in qi may result in restrictions on c 
1 


description as c' Since it may be necessary to change 651 from the transfer 


In any event, c,, assumes a new 


11 1 


11° 


function dictated by non-interaction theory, it now becomes 651- Therefore, 


፲..6፡.09ሠ፡0 ባቸ (4-15) 


The procedure is to make 1 t 91 1 8° that there are no zeroes in the right 
half p plane. Note that 1 + dii will have zeroes in the right half p plane because 
dı = “ሃራ can not be properly compensated so that Z = 0. 

In as much as the aij are fixed by the physical system it is much easier 
to design e and/or 1 So that there are no right half p plane zeroes in the 
resulting 1+ di AS a first approximation, write 65) in the form dictated by 


non-interaction theory, which is: 


Ca = — r Ca (4-16) 


22 


In as much as c! 


11 is arbitrary in some sense, make it as simple as desired. 


Then substitute into equation (4-15) and apply Routh's criteria to the resulting 
polynomial to determine whether or not these are zeroes in the right half p plane. 
If this test fails, alter 651 by a factor, a, where a is a polynomial in the operator 


p or a ratio of polynomials. Again, test the resulting polynomial for positive 
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roots. Eventually, a q will be chosen such that 1 + di has no zeroes in the 
right half p plane. In addition, if q is chosen so that the "break'' frequencies are 


above the system natural frequencies then the interaction will be restricted to 


these higher frequencies. | 

Having removed the difficulty in dip attention can be focused on doo: An 
attempt should b& made to design it in accordance with the general theory of non- 
interacting systems. If no compensation difficulties are experienced then Coo 
and Cio will be dictated by the general theory. In addition, 1+ d, -1+ bo 
and there will be no right half p plane zeroes. If compensation difficulties are 
encountered then follow the procedure outlined for d p 


Naturally, the system has some interaction which appears through the 


915 and d; 1 transfer functions. In fact, for the case presented here 


€ ን 
di; = Qi er + Aiz 2 = © [E a 2 (4-17) 


da) = Gi + QC, #0 (4-18) 


Due to interaction, the overall system closed loop stability is dictated 


by the zeroes of: 
۱ / / 
(Led) JO diz) - di AO (4-19) 


Determining the stability through this equation has already been discussed at 
length in the previous section. In the event that there were zeroes of equation 
(4-19) in the right half p plane they would stem from the interaction terms, 


dj» and doy: Thus, it would be necessary to choose another q factor so that 
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both the 1+ d and I= dz dh: contributed no right half p plane zeroes to 
(1+dn Vi tdia) 
equation (4-19). 

In the case considered here da= O. Thus, system stability is assured 
if )۱+ di ) ( I+ de have no zeroes in the right half p plane. Actually, 
a = 0 due to the slight deviations in 2 and E of equation (4-17). How- 
ever, it is not expected that such slight interaction through this term would cause 
instability when coupled with ds, 

In summary, when compensation difficulties arise and interchanging the 
rows of the ፅ vector does not alleviate the problem, then it is necessary to 
deviate from the non-interaction theory to the extent that there are no zeroes 
in the right half p plane of the resulting characteristic equation of the system, 
namely: 


(1 td, )((፦ጩ) As de O) 


| 
Naturally, interaction will be introduced through the transfer functions. In 


the case considered here, the amount of interaction as a function of frequency 


can be found from the following transfer functions: 


4x, = 3 (4-20) 


dz, (4-21) 


f/x, DT 


D| = (i +d IG *dii) adi. 4, 


| 


where 
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Chapter V 
The Synthesis of an Aircraft Control System 

In present day aircraft, the auto pilot performs many functions in the 
integrated control system. Although the system inputs are generally restricted 
to combinations of elevator, rudder, aileron and throttle, the system controlled 
outputs are many? For instance, an auto pilot system could control the aircraft 
velocity vector through the control of the pitch, yaw anb roll rates. Other 
systems could control either altitude, angle of attack, airspeed, glide path, 
Mach number, lateral acceleration, turning rate, etc., or combinations of the 
foregoing. It is the intent here only to point out that the auto pilot can be used 
to control many variables and each subsystem of the auto pilot can control as 
many outputs as it has inputs. 

In order to point up the theoretical analysis and the practical problems 
associated with a design based on the general theory of non-interacting controls, 
an auto pilot: subsystem shall be synthesized in this section. Naturally, the 
goal of such a study is to arrive at a controller design that permits good system 
performance as well as the desired non-interaction. However, as was pointed 
out in the preceding sections, a perfect non-interaction over all frequencies is 
not obtainable nor is it necessarily required. In addition, complex physical 
systems such as an aircraft could easily require impractical or impossible 
controller compensation when applying the general non-interaction theory. 
Thus, the designer may have to accept some interaction in order to preserve 
system stability. Furthermore, he may be forced, at least in part, to deviate 


from the general theory in order to design a practically realizable system with 








31 


desirable characteristics. 

The subsystem to be synthesized is one dealing with longitudinal control 
in which the inputs are elevator and throttle deflections. The controlled out- 
puts are glide path angle, 6, , and indicated airspeed, /Up . It should be 
emphasized that this study was chosen only as a means of lending practical 
significance to the theory and is not considered as a demonstration of a specific 
auto pilot design. 

The assumed aircraft and flight condition is an F-86 '"Saberjet' making 
an approach down a 3? glide slope. Small perturbations from symmetric 
equilibrium flight will be assumed thus permitting the use of the linearized 
longitudinal equations in the ensuing analysis. Upon determination of the 
controller design for the best system performance, the synthesized system 
will be studied on the differential analyzer. Such a study should clearly demon- 
strates the applicability of the theory, error analysis and the practical system 
problems. 


The large motion longitudinal equations of motion for an aircraft using 


a combined body- wind axis system are shown in Ref. 5 to be: 


Z Fy =0 ==mWp Fw 9 —mVp tl cosa + Pasino (5-1) 
Ex AS 


ZFwa = 0 = MVpQu +Fya = ۷۸۷۵ 40۷ - 2 5۱0۵۷ + Ba cosol (5-2) 
+4 mg cos Ow 


(973) 


2 My -0 =-244 9 224-12 ۳ LI [REP ተ ኮጎ 
+ Ty + Toq (P+QR) +Iy2 (R=PQ) 
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The linearization of these equations is carried out in Appendix I. In addition, 
a fourth linearized equation is obtained so that there are four simultaneous 
equations involving the four output variables, Ne : 5 ; Y : 6 and the inputs, 
8 and Du | 
It is desired to determine the transfer functions that relate Np and 7 to 
2 and x . It should be clear from the theory that these transfer functions 
are the ai elements of the À (aircraft) matrix. The mathematical steps taken 
to determine these elements from the four simultaneous equations are completely 
outlined in Appendix I. 


The elements, a;;, are merely a ratio of polynomials in the operator p 


ij 
or D. However, they derive their existence from the coefficients of the line- 
arized longitudinal equations. These coefficients are composed of thrust terms, 
gravity terms and aerodynamic force and moment terms expressed through the 
aircraft stability derivatives and trim lift and drag coefficients. 

The stability derivatives for the F-86 were obtained from Ref. 6 for 
flight at sea level at 163 knots. They are listed in Appendix II along with 
other physical data pertaining to the F-86. Calculation of the trim lift and 
drag coefficients as well as the b; and i coefficients of equations (1-9) 
through (1-12) are also found in Appendix II. Knowing the bij and fij co- 


efficients permitted the calculation of the ai elements from equations (1-18) 


through (1-21) of Appendix I. The results are: 


a "رر سے‎ (D-.012N D*+ 1.422D +3.4771) 
መ... መሠ 


(D*°+1.4144D + 3.4534) ( D*+.01S0D +.0245) 





33 


(D* +14 144D + 3.4534 )(D*+. 0156 D +.0245 


ETT (D-150)(D+1.1) 5 
aa ) 


77 c ۳۳ D+ 1.409) (D*+. 701 D +3.087)_ 
) +1.4144D + 3.4534 (D+0156D +.0245) 


= 07 (D-.00717)(D+7.05)(D-8.5) 3 
(D +1.4144 D +3.4534)(D*+.0156D +.0245) 


% c ERR 


The two quadratic factors in the denominator of the a terms provide 
information as to the period and damping of the short period and phugoid motions. 
The period and damping factors obtained are: 


Short period motion Phugoid motion 


EU xm ہے‎ RI um a A 
iB m = 3.38 sec. pua 40 sec. 
f. سے‎ -582 tp = . 0842 
It may be seen from Figs. 1 and 5 that in vector-matrix notations: 
y - 55 where Î = (tii) 
NE CS) 
A - (au) 
For the specific problem under study: 
Yi = Vp ; 4 7 "e ) ۷۳۹ dr ) ES 
o En ^ o C c 0 
ولا‎ > Ow 2 بت ولا‎ ) à = de ó = de 
Therefore: 


Up = (ti Qu Su) & + (tn On Saz) de = Gy dp + Che de 


(5-4) 


On = ( د‎ Az S +l(t 3 $22 ) de _ Qe, d+ + de 


(5-9) 


To determine the aij terms of the above equations, itis necessary to specify 
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the elements of the diagonal servo and transducer matrices, S and T respectively. 
This has been done in Appendix II. In addition, reasons governing the specifi- 
cation of these elements or transfer functions are included in the appendix. With 


this knowledge, the aij terms of equations (5-4) and (5-5) are found to be: 


(DP OED Ha S4771)‏ ےہ 
(D* +1.4144D +3.4534 ) (D> +.0156D +.0245) (2D+1) (.2D+1)‏ 
ıa = _,gos| _ D-1So) (Dt. I E‏ 
Dt +1.4144D +3.4534)(D'+.0156D +0245) (04DH) (2 D+1)‏ 
ie, = .00625 D+1.409) (D*+.701D + 3.087 -‏ 
.=| 
(D*+1.4144D +3.4534 |( D'0156D +.0245) (aD+1)‏ 


(D e1.4144 D € 3.4534) (Di. 0156 D 0245) (.04D+1) 
It has been pointed out that d. z 6 /8.. is an open loop transfer function 
ii 


for one "channel" of the auto pilot system. From Fig. 5, it may be seen that 


the open loop vector-matrix equation is Å = ACE. It follows that: 


d, = Yo, = ታሃሬ = lag! = Cy JAL = PE. (5-6) 
7 vr 
dzz = Vo, = de/e, = Cze 0 = a = 2 (5-7) 


It is evident that equations (5-6) and (5-7) are the open loop transfer functions 
for the velocity and pitch angle "channels". It is desired that these channels 
be closed loop stable as well as having desirable response characteristics. An 
attempt is made to design the controller transfer functions Sr and oo so that 
these requirements are met. If they are satisfied, then the other controller 


transfer functions, 515 and Cy, are determined from equation (2-14). 





Controller Aircrat# Trans ducer y 


@ | A 7 


Pie. 5 The Open Loup Control System Configuration 
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The determinant of the A matrix was found to be: 


ES 


| a, 01۱ | 


ni D+7.74)(D-8.587)(D* +1. 41S D +3.456)(D'+.0187D +0214) 
( 


|۸ 
2D) 2D) 04D) (D^-1.4144D 3.4534) (D 4.0156 D 4.0245 


Substituting into equations (5-6) and (5-7) and cancelling equal or 
essentially equal terms in numerator and denominator (excepting terms in 
کله‎ 


right half plane ) provided: 


7 0 ]لہ = % 


00117 p-8,5) (2 p+!) (-2 pH! ptt) (8-8) 


7 (ا-م65(.) (ا+صؤد).)  ا‎ 
lo, = lag (هده.)‎ e 238 p"+. 409p+1)(04p+ 1) E 


du 


In deciding what c, Should be in order that 
+ dy 


be closed loop stable 





with good response characteristics, Nyquist plots had to be used. A sketch 


Yo 
۹۹3 


half p plane of ا‎ > Nyquist stability criteria indicates that N = + 2 for 
H 


of the Nyquist plot of —- is shown in Fig. 6. There are two poles in the right 


Z-0. This implies that c, , must be chosen such that the Nyquist plot of the 


11 
resulting © or dj, must encircle the critical point, -1, 0, twice in a counter 
clockwise direction. Without doubt, the required compensation presents a 
rather formidable problem. In fact no practical compensation E be found. 
Thus, itis clear that the system will be closed loop unstable in that it can not 
be properly compensated by 11- 

In an attempt to alleviate the above design problem, the rows of the 7 


vector were interchanged. As a result, it was possible to compensate dii in 


accordance with the general non-interaction theory. However, compensation 
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problems arose in o or doo and no practical means of compensating were 
found. Thus, this method failed. 

The next approach was to alter 91) with a factor q such that there were no 
zeroes of lt dt 1 in the right half p plane. This approach is discussed at length 
in the section on system problems. From that section, it is known that a more 
general form ell + x or 1 + d is given by equation (4-15) which is repeated 


l 
here: 


ltd, = Qu OE M (4-15) 


It was pointed out in Chapter IV that c! 


11 9T 651 could be somewhat 


arbitrary when di 1 could not be properly compensated by the non- interaction 
design methods. In the case at hand, instead of making 511 completely arbi- 
trary, let us assume that the p-8. 5 pole is cancelled by the p-8.587 zero in 
equation (5-8). Of course, it is understood that this is incorrect in that one 

of the right half p plane poles are being neglected. Nevertheless, compensating 


the remaining o with P = 1 results in 51 , being: 


a ۳ 2 epti 
5 Sa "i (5-10) 


Naturally, this is not the necessary compensation for closed loop stability of 
du 
H-du 
AS a first approximation to 651) let it be dictated by the equation: 


as dictated by non-interaction theory. 





Ce SS ረረ. e 


22 


Substituti í in: 
bstituting for 291» 255 and C1 results in: 
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C/ = 00292 (+«711pt1)(.04 p+1) (324 på +.227p + I 


B. (^ !p*i) Cp 7.00717) (430€ pr 1) (^ 118p - 1) (o 


Substituting equations (5-10) and (5-11) as well as aj, and ajo into equation 

(4-15) and applying Routh's stability criteria to the resulting ninth order poly- 
له‎ 

nomial indicates that there are zeroes in the right half p plane. However, by 

altering دوه‎ above by a factor do» where qy = -(. 1178p - LY 04p + 1), the 

new Coy becomes: 


2 
e D 0022 (TIPPET DE 227p+1) (5-12) 


(p-:00T17)(.1p+1) («1304 p+1) 


Going through the same substitution procedure and applying Roth's criteria 
indicates that there are no zeroes of 1 +dy1 in the right half p plane. It follows 


that the velocity "channel," Y,, is now closed loop stable. In addition, inter- 


1 


action is present in the system because the 511 and %1 of equations (5-10) and 


(5-12) are different from the c}; and c3} that would have been dictated by 
following the non-interaction theory. 

The doo transfer function was determined from the non-interaction theory. 
It is shown in equation (5-9). A Nyquist plot of رہد‎ is shown in Fig. 7. 
By making Coo eu su ; d, o (or Y o5) was properly compensated so that 
there were no zeroes in the right half p plane. Furthermore, the compensation 
was chosen so as to give good transient response. 


Having determined c. in accordance with the non-interaction theory, then 


22 


it was permissible to find Clo from equation (2-14). Therefore: 
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Co = al - - 2.045 (1:2p+1)(.Ilpr !)(zprı 7 -1) (5-13) 
Ay (p--012)(o4p+1)(o2p+1)( 288p*+.404p+1) 


In as much as terms such as (. 00667p - 1), (. 02p 4 1) and (. 04p * 1) are 


only effective at the higher frequencies, c, , can be simplified by neglecting 


12 
them. Naturally this will give rise to interaction at higher frequencies. But 
this is of no real consequence provided that these gross deviations do not im- 


pair the system stability through the im 1 term of the system characteristic 


2 
equation. In the problem considered here the stability is not impaired by the 
interaction terms. 

In light of the foregoing: 


/ | )/ خم 91.)( لومت‎ 
m = 4 C = 905045 (^ept CUR er!) E (5-14) 


(P--0l2) 288 p*+, 409p +1 ) 
as = (:04p+1) (02 p+1) 
| - (0067p — ı) 


Attenuation plots of e and the c,. of equations (5-10) through (5-14) 
1) 


“22 
are shown in Fig. 8. One may readily see that the effect of neglecting the high 
frequency terms is to level off the gain at high frequencies, such frequencies 
being above the natural frequency of the overall system. 

The amount of interaction contributed through d 9 and doy can be found 
from equations (4-20) and (4-21), In the case considered here, do and u. 


are as follows: 


e f GE + 2 Coz (5-15)‏ عله 
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/ / / 
de, = تت‎ Cer Uzz (5-16) 


where c!,, is given by equation (5-12) and c is given by equation (5-14). 


2l 


Computer Simulation 


12% 


With the 51, designed, the auto pilot system was mechanized for study on 
the Reeves A Analog Computer (REAC). The mechanization used is 
shown in Fig. 9. 

System response to step inputs of command velocity, Jp: , and command 
glide path angle, Ow , are shown in Figs. 10 and 11. Note that the interaction 
is essentially zero except at high frequency. However, this was expected in 
that high frequency terms in Er and E were neglected. There is, although 
not evident in Figs. 10 and 11, a low amplitude damped oscillation that affects 
the steady state character of 5 and EN . These oscillations are the direct 
result of disturbance inputs, Uy , that enter the system through the E In 
fact, Figs. 12(a) and 12(b) show these oscillations in ግ O, when there 
are no inputs other than disturbances in the system. 

As a means of verifying the magnitude, frequency and damping of the low 
amplitude oscillation, a step input, Xo, Was introduced into 521 through the in- 
put to amplifier R-8 of Fig. 9. The effect of the disturbance input on p 
6y is shown in Figs. 13(a) and 13(b). 

An analytical study of the system was made in order to determine the 
relation between E and x, and also to verify the traces seen in Fig. 13 as 


2 


those that stem from disturbance inputs. In the case where there are no desired 
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perturbations from steady state condition, then م لہ‎ =0. Therefore, it may 


be shown that: 


Az Ca 
Nor پر‎ 3 
P | + Au Gy + Az Ca : e 


Substituting into equation (5-17) and simplifying results in: 


Npa (24 DG7U pt) (pt tet) (p-1S0)(324p"+.227p +1) - 
P [p'+ 44.۱ p® +6۵۵ *م 5۹70 + 7م‎ + 27 0SDp° + 90,400 p* + 131600 p> + 71,000 p*- 
+5150p +4100] 


(5-18) 


where the static gain is unity. 
Application of Routh's criteria to the ninth order polynomial in the denomi- 
nator of equation (5-18) indicates that there are no positive roots. Therefore, 


the system is stable for the disturbance input, x,. In addition, it may be shown 


2 
that + - 0194p + . 0210) is a factor of the ninth order polynomial. Further- 
more, it is the only factor which has a low characteristic frequency. In fact 
this frequency is .0231 cycles per second and the damping factor is = .067. 
It is quite apparent that this quadratic factor alone accounts for the lightly 
damped high period oscillation seen in Fig. 13. Naturally, similar analysis 
can be made for other disturbance inputs into the system. 

The low amplitude oscillation discussed above as well as other responses 
to disturbance inputs can have a detrimental effect on the system due primarily 
to their magnitude and damping. The addition of a high gain loop would introduce 


gain prior to the introduction of noise and hence make the system less sensitive 


to the disturbance inputs. 
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Chapter VI 

Conclusions 
1 In theory, control systems based on non-interaction of control variables 
are feasible but somewhat complex. For complete non-interaction, the theory 
requires precise mathematical knowledge of the behavior of the system to be 
controlled and exact mechanization of the controller transfer functions. In 
reality, these requirements generally can not be met but can be approached 
to the point that the resulting interaction is negligible for all frequencies. 


2. Small deviations, €; , in a,, due to insufficient knowledge of the physical 


L J 


system and/or in ei due to the inability to mechanize transfer functions exactly 
should not impair the stability of the system if the dj or open loop transfer 
functions, are designed with a reasonable amount of stability. 

3. Gross deviations in E so as to obtain a simpler mechanization can 
result in considerable interaction and possible instability in the closed loop 
system. If deviations are restricted to terms effective only at frequencies 
well above the natural frequency of the system, then the interaction will be 
restricted to these frequencies. However, system stability can be impaired 
through the interaction terms now present in the characteristic equation of the 
overall system. This equation should be investigated for zeroes in the right 
half complex plane. 

4, System errors due to disturbances entering the system can be large. 
However, these large errors can be reduced considerably by addition of a 
suitable high gain loop. 


5. The system error equation indicates that the overall system can be 
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stable for command inputs, X; and yet be unstable for the disturbance inputs, 
Uy and Ne . In any event, a complete error analysis for disturbance inputs 
should be made. 

6. Zeroes in the right half p plane of the co-factors of 44 can create compen- 
sation problems through the open loop relation, اه‎ = Cü . Naturally, €; must be 
chosen so that there are no zeroes in the right half p E of 1+ 91. Making 
the proper choice of : can be a formidable problem in some instances with 

the result that no practical تحت‎ can be found. 

ds Interchanging the rows of the control vector, § , will alleviate the 
impractical compensation problems in some cases. 

S When compensation problems are encountered and changing the rows of 
the ፅ vector does not alleviate the problem, the non-interaction theory can 
not be used for the design of that particular "channel" or open loop part of the 
system which is causing the problem. Some other approach must be used to 
make the "channel" closed loop stable. Normally, the end result of deviation 


from the non-interaction theory is system interaction, which, it is hoped, can 


be made small. 
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EU PENDIX T 


LINEARIZED LONGITUDINA 


L EQUATIO 
OF THE Elm 


NS AND DETERMINATION 
NENTS OF THE 


A MATRIX 
A. Linearized Longitudinal Equations 


The equations for longitudinal motion in the 


general 
or large 


motion case ina combined body-wind axis system 
have been shown in Ref, 5 to be: 


(I-4) ae O = -m Vp + Fx z—mVp4-R casa +f Sind 

TX = DE Ou 

ይ 2) E 'we 207 «mV Q,, +Fwz = mVeQw — E sinat - P cased 
+25 + mg COS Ow 

۷ 7 Te ፌሪ 

(I-3) 2 My =0 =Lyy @ + [T+ JER + Iex [R= P] +- 

Considering small perturbations from a steady state 


longitudinal flight condition, the above equations can be 
linearized as follows: 


1. Assume the following 


COnditions: 
P-R = بے‎ - ٠. 5 Lye = 4፡7 = 0 
O= Oy = - 9 O 

2 


. Assume the following steady state plus perturbation 
relations for the system variables: 


Vp = ኣቱ + % 
و دوک ات‎ f 
Ow = Ov, + 6, ን Or 
Q= Q, ፦ሩ 
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where Up, a, Bu 9, Gu Pe AT are the perturbation 
variables. 


5. The dimensionless coefficients, X,,25,Y,,N,, lg, and 
Mg are functions of xx ©, fs, Qs and Rg» By using 


non-dimensional angular rates, X P, Q and pM 


then X Zg and Mg become functions of X, ኃረ a 


g? 
anā control inputs only. Thus, neglecting the lesser 


inportant derivitives, the dimensionless coefficients 


become: 


(1-4) X, = EVS | Gx +C] 
(1-5) Zs = po s[Ca +Ca +Ca 7 « Ca, Se] 
(I-G) Ms - L o VpS | Cu Cuz e ^ Cy, + Cm Mz Q 
| + Cm, E 
4. In the wind axis system, it is true that: 
(r-) Eu = Qu 
(r8) & = Q-Qy 
5. Linearizing equations (I-I) through (r-3) about 
the steady state condition (subscript one) and 
rearranging in operator form with D=d/dt provides: 
(1-9) b, wp + ba X +b, Ow tba q =f bp + Se de 


Where : 
D 


| »6 می 5 9 4 0 














z= سح‎ «V = | E ane = E 5م‎ Cx, 


J m di. 
bs - 1 = NV 


i 





| 


Fu 
Je 


Ó 


۱ سا وت 
m: (255) COS, Ho sn |‏ 


ያ‏ : ہے ہی 


( Z-ı0) Da Up +|; << a2 Dio 2 , tb 1 = VER = + E. 3 
uo ! 


ba = Koy = 


"zv 


ጠ 


I ይ + Cao Se | 

be = D= Kay = D- ۱ | ری گام‎ Beast, Bom] 
:د‎ ቬ | sin Ow, | 

=~ [ike] = - [1+ ESE Cag | 

ae‏ ا 


(I-! 





l) boi + بط‎ bb ر کے مو‎ 
whae: 
1 ons su 
ba, = - Kqv ን መመ | SES [Ma ©; + [mo Hm; fe | 


Vo, >‏ ج نا 
DE =ፎ Bio en a‏ سے سے 
e|‏ : مه b32 1 D EK |n] 25, ር? t‏ 


DN NT 


سب ما 


ba = D-Kgg «[P- هه‎ Ca | 


ዛመ--- 


p = 


O Vo, JC 
پل‎ = Kase = E Cms, | 
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The fourth equation is found from equation (1-7) 


aná (I-8) as follows: 


Shy = (Qu od which when linearized 
with Ow, zd = Q, E provides: 
ፎኔ + ዉ ፦ 4 =0 


Fearranging and using operational form results in: 


(I-12) ba Ap +bae 5 +bas Õu + baa g = fn Srt fn de 
hore | | 


ba ہے‎ 

542 = 13 2 +/ك‎ 
D43 =D zd/ge 
ኤፊ = ፦ | 

ተሩ) = ርጋ 

28 = O 


B, Determination of the Elements of the Aircraft Matrix 
The vector-matrix equation for controlling airspeed, 


መሙ 


/ , and the glide path angle, Qy, using power and elevator 
inputs, 2 and 7 1s: 


y ros 0 5 = ۱ [ E 
ا‎ 6 - ٣ ۳ or "OD 
ENG 1. (s m 
When the relation above is expanded, the result is: 


(I-14) op = ay, 1 t ፈሬ ðe 


(EP Su -.። 7 a T em A 


The elements, @ij , in equations (I-I4) and (r-!5) are 


found by solving equations (I-49) through (I-:2) simult- 
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aneously for Vp and Ó, . Using determinant methods yields: 


M b | Ö bi biz $F, E €. 
— (f, = + Fade ) på baa D24 — ba; bii a 7 ነፈ de) bza 
ملم‎ E de bar O ba Ow = [ba جوا‎ ፌሪ. ودد‎ 
0 baz عمط ۵ الذي دوم‎ 5 E 
A A 
where | bi bag | . O 


Da 9 ፓም‏ = ھ2 
ba ወጩ bay‏ روط 
ات O baz bys‏ 


Expanding and simplifying the notation results in: 


5-%) %» BWA-A OJS , 7 


bs لیوط‎ baz = ba baa bi | 


N A 
(T-17) o [&(E) - &(F)]& | [&,(6) - &F) | e 
A A 
kur: = 
A) = | bzs bsa baz tbri bas bas + Daa bas — bas bas baz | 
= baa baa + baz — bag bea bi | 
)٤( = | ba4 bea + baa - ba bas «bys bas be] 
(E) = 
) 


——- 
سا 


| 
E 224 ba b42 Tb bs = ba ہے‎ ba; b34 2 
| 


å De — y را‎ — bar P T 
A =| b, (A)- ba (B) Gee (C )] 
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From equations (1-16) and (1-17) it is clear that the 


matrix elements are: 


(I-18) 


(I-19) 


(1-20) * 


(I-21) 





Bo. áA- fa, (0) 
ii 290 
Qa, = ጋ (C) fr | (55) 
A 
O = ROE A (EF) 
aNd 


ts, (G) E hi (F)‏ = بل 
AU‏ 
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APPENDIX II 


CALCULATION OF THE F-86 STABILITY DERIVITIVES AND 
DETZRMINATION OF THE COEFFICIENTS, bij AND fa 4,0F APPENDIX I 


A. Physical characteristics of the F-&6. 


The following data has been obtained from Ref. © . 


S = 288 ft? 


W = 15,500 lbs. (lightly loaded) 
b = 357.1 ft. 
6 = 8.09 ft. 


= 8625 slug-ft? 
XX 


2 
yy 00 slug 


x: 2 


بر ہم ہر ہم 
} 


slug-ft?‏ 35300 = ر„ 
B. Assumed flight condition.‏ 
Altitude = sea level‏ 
Mach No. = .246‏ 
۳ / و 0055258 = density‏ 


= 163 knots = 275 ft/sec 


Q-4- 51 = OX, = Gu, EO‏ ح ولا دول 


0 7 
EN = — 5 = —, 0523 radians 


The above data is for an F-86 proceeding down a three 
degree glide slope at 163 knots in a "clean" configuration. 
C. Lift and Drag Coefficients at Trim; Stability Derivitives 
RW COSOWw 


7۶7 اسما 
Cups‏ 


l. Ce, = 
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2. From the drag polar using Cr = Si7as the entering 

N 

argument: i 
5 B 


= .04-45 
Therefore, ከ ከ ከህ = 1152 1b3, 


3. Stability Derivitives 


The following stability derivitives were obtained 


directly or indirectly from Ref. © 


ገ 
Cg ,-.065 per degree = 3.72 per radian 


C - مه‎ ace = ¿J01 er radian 


7 Sch 2 — ~.433 per radian 
ma = SE SE 
ር) =-5.6 per radian 


2C 
(Da 2 26; E = .51 per radian 


per radian‏ 1 > و ما 
Cho -1.08 per radian‏ 
por radian‏ 495 = وت 
Cy. = -+1775 per radian‏ 
oa‏ 
Cu.= CUms. = -.95 per radien‏ 
ia. = de‏ 


E و اند‎ SCE PEF ean 
ጆ 3. ው CL Se va 
fficients of longitudinal equations 


The addinonal assumptions are | 
de, I E = 7, =10 

5 Vi 5 ^ ^ 2 ስ 

0 Ky ኑ CP کس = و‎ 4 


trim = — 00054 ALL 


۷0 WR S ft 





E. 


ዙ‏ م 


3 


K Vol 


K ved = 


(Kus, 
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- pes 
EL ረፌ k- = Å DDE ፈጋ 
—. 6067 per ad an. 
= d 
۳" E 
= Paes cos. راح‎ 
ing (Bip) CO 318 
= 25 |, 3 
ا‎ Ca] = ESEC, n] = onos £ 


dom 
Vp L ie S Ca,- Fx cosa څ‎ 


IEL = j- ፡ 





E E 
L 
= E 5۱۲6 = +. 00613 بط‎ 
= 7 
E E| Ca] = E. 2፡1. د‎ 
= rad. 
- E28 s 
uS Y udi Z 
| 2 i E ے٣٦‎ 
= د ,ملام‎ Ca I 58 / 
ر لل‎ de ሮም” ር ኃሬ ብበ 
= Oo ? ری مهه مر‎ 1 pns 
= pvp Se” Em y = —,//3 ل‎ 
۷ 
B مام‎ 7 => c dot 
> 2) 
= ሪብረ ربج‎ = 58? — 
27 N are 
7 
ws 
erp: = Cos, = — 6,74 —_ 





a 





D. Coefficients of Equations (I-9) through (I-12). 


b,, = D/e - Ky = +0311D + .00054 


0 


سن رب 
pop‏ 
مر سم 
ዘ" |‏ 
ዘ‏ 
® 
E‏ 
Fr?‏ 
CO‏ 


Kos 
boj =-K,y = .00085 


boo = D- K, =D + .709 


Dog o, =-.00613 


daa =-[1+K,,] =-.9853 


የ21 = ፻25. = +.00624 


r | 
Dya = = |DKq + Kog | = ۰1130 + 9 


ከ33 = 0 

byg = ከ. = DES 
pO 

Ras, = 4‏ = وی 
Det‏ 

=D = ٤‏ ويه 

Das = z d 

٠‏ دو 

لا 1ے 

PE. 
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APPENDIX III 
ELEMENTS OF THE SERVO AND TRANSDUCER MATRICES, T and S 


Å. Servo Matrix Ele ments, s, 41* 


l. The relation between elevator deflection, A , and the 


signal to the elevator servo, ሪኃ , is defined as 854° 


A typical elevator servo system can be described as a 
second order system with natural frequency, Wh, , and damping 


-ኃ 
ratio, te y The phase relation for such a system is given 


N 


by: (2 Z leùn 
۱ ME "ne 
OI) Ø =-tan" I ከ. mig 


Ne 
UT n the closed n transfer function has a db 





versus log w characteristic as follows: 


magnitude 


(db) 





It follows that if the elevator servo is designed such that 
the natural frequency, Un, ‚ 15 larger than the expected 
range of system frequencies, Ay, or W<< Un, then equation 
(III-1) can be approximated by 
۲ ۱ 1 2 

(111-2) Ø =-tan = ٧ 
If the damping ratio is taken as .S,then the phase relation 
becomes: 

IIIS) 2 = - ተፅ.” Y = dun ED 


۱ 
where Va — 
Wing 





It may be shown that the phase relation of (III-3) is 


that phase relation for å simple first order system of the 


| 
SB ST 
(III-4) (p? Ir Tp! 
n 


In conclusion, if the damping ratio is chosen as .5, 


form: 





1 


and the natural frequency of the servo is taken as 25 rad/sec, 
Y 

then"the expected range of system frequencies the second 

order elevator servo system may be approximated by the first 


order system: 


De |‏ 
— — ی تج ار ms‏ 
| ان 4 . "m‏ ; 


2. The relation between effective throttie perturbation or 


thrust, 5, , and the actual throttle perturbation, dr, is 


simply described by the following first order system: 


7 | 


(III-6) تچ سس‎ — 
Ör Tor 


where ‘T= the time delay in obtaining the new 
engine RPM which is representative Qf the 
thrust. 


The above relation is pointed out in Ref. 7 and is 
generally true for most turbojets. In addition, a typical 
time delay of two seconds, T=2.0 , was obtained from this 
reference, Thus: 


(III-") m 
SPE or | 


mie SE Sy SE مسح ع‎ | 


ፋ 2p+ I 


B. Transducer Matrix Elements, 5414 * 


l. The transducer that relates actual indicated aircraft 





velocity, E to the measured aircraft indicated velocity 


is given by: 


(፲፲፲-6) t,(p)= E | 
^p eps 


በ 





Such a transfer function assumes no position errors in the 
airspeed indicator nor any dynamic errors in the instrument. 
The time delay, T= .2 sec , is due to the pressure trans- 
mission delay in the pitot system. 
The basis of the foregoing statements is derived from 
a discussion of airspeed indicators and measuring devices 
contained in Ref. Å. 
2. The transducer that relates actual glide path angie, 
e to measured glide path angle, On , is defined as tone 
It is assumed that the angle © and the angle of attack, o 
are obtained without error from a gyro and an angle of 
attack indicator. These two quantities are added together 


to produce وت‎ ۱ 5 


)111-9( 





eus (ps ——-4 
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